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Wavelength  dependence  of  laser  ablation  of  silicon  was  investigated  with  nanosecond  ultraviolet, 
visible,  and  infrared  laser  pulses  in  the  irradiance  range  from  3  X  10ll)  to  1  X  1012  W/cm2.  For  266 
and  532  nm  laser  pulses,  the  depth  of  laser-produced  crater  shows  a  dramatic  increase  at  a  laser 
irradiance  threshold  of  approximately  2X1010  and  4X1011  W/cm2  respectively,  above  which, 
large  micron-sized  particulates  were  observed  to  eject  from  the  target  about  300-400  ns  after  the 
laser  pulse.  In  contrast,  for  1064  nm  pulse,  this  dramatic  increase  was  not  observed.  The  underlying 
mechanism  for  the  observed  threshold  phenomenon  is  presented  in  this  study,  which  can  be 
attributed  to  the  thermal  diffusion  and  subsequent  explosive  boiling  after  the  completion  of  the 
interaction  between  the  nanosecond  laser  pulse  and  silicon.  Based  on  our  delayed  phase  explosive 
model,  the  ablation  depths  were  calculated  for  different  wavelengths  and  compared  to  experimental 
results.  Plasma  shielding  during  laser  irradiation  was  included  in  the  model,  which  plays  a  key  role 
to  the  coupling  of  laser  energy  to  the  irradiated  material.  ©  2008  American  Institute  of  Physics. 
[DOl:  10.1063/1.2978369] 


I.  INTRODUCTION 

During  the  past  few  years,  interactions  of  laser  pulses 
and  materials  have  attracted  considerable  attention  due  to 
their  applications  in  a  number  of  areas,  such  as  deposition  of 
thin  films  and  laser-based  chemical  analysis.  ’  One  key  pa¬ 
rameter  in  characterizing  how  well  is  the  coupling  of  laser 
energy  with  a  solid  is  the  relationship  between  the  ablation 
depth  versus  the  laser  irradiance.  The  removal  of  mass  from 
the  target  by  a  laser  pulse  can  occur  by  both  thermal  and 
nonthermal  mechanisms.  For  nonthermal  mechanism,  it  is 
suggested  that  the  incident  laser  radiation  can  induce  a  large 
population  of  electrons  to  a  highly  excited  nonequilibrium 
state  near  the  surface,  which  leads  to  bond  breaking  of  the 
target  material  and  subsequently  causes  ions  ejected  from  the 
target  surface.  For  the  thermal  mechanism,  the  excited  elec¬ 
trons  transfer  absorbed  energy  to  ions  by  the  electron-phonon 
relaxation  process,  while  the  heat  is  conducted  into  the  solid 
lattice.  The  conducted  heat  can  melt  the  sample  and  bring  its 
local  temperature  far  beyond  the  boiling  temperature.4  For 
nanosecond  laser  interaction  with  solids,  thermal  mechanism 
appears  to  be  the  dominant  mechanism.5 

When  the  laser  irradiance  exceeds  109  W/cm2  for  a 
nanosecond  laser  pulse,  the  temperature  near  the  surface  of  a 
solid  target  can  exceed  the  boiling  point,  thus  a  superheated 
liquid  layer  is  formed.6  It  was  suggested  that  explosive  boil¬ 
ing  takes  place  where  homogeneous  vapor  bubble  nucleation 
occurs  when  the  target  material  reaches  ~0.90Tlc  (Ttc  is  the 
thermodynamic  critical  temperature).  As  a  consequence,  the 
target  material  makes  an  abrupt  transformation  from  super¬ 
heated  liquid  into  a  mixture  of  liquid  droplets  and  vapor. 
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which  are  then  ejected  from  the  target.  ’  Laser  ablation  from 
single-crystal  silicon  with  laser  irradiance  109— 101 1  W/cm2 
(single  pulse)  was  studied9,10  in  an  experiment  using  a 
Nd:YAG  (yttrium  aluminum  garnet)  laser  with  266  nm 
wavelength  and  3  ns  pulse  duration.  The  results  showed  that 
the  ablation  depth  increased  dramatically  at  the  laser  irradi¬ 
ance  threshold  of  approximately  2  X  1010  W/cm2;  above  this 
threshold,  large  particulates  were  ejected  from  the  target  sur¬ 
face  about  300-400  ns  after  laser  interaction.  A  laser-induced 
transparent  layer  when  the  temperature  approached  the  criti¬ 
cal  temperature  was  speculated  to  result  in  the  abrupt  in¬ 
crease  in  the  ablation  depth  at  the  threshold.10  However,  ex¬ 
periments  were  unable  to  verify  such  a  transparent  layer 
during  pulsed  laser  ablation  of  solids. 

We  developed  an  alternative  theory  based  on  delayed 
phase  explosion  model"  that  was  able  to  interpret  the  experi¬ 
mental  observation  of  UV  laser  ablation.  In  this  article,  we 
extend  our  theoretical  analysis  by  examining  laser  wave¬ 
length  dependence  of  experimentally  observed  threshold  be¬ 
havior  of  phase  explosion.  In  the  following,  we  first  present 
an  overview  of  experimental  results  of  ablation  depth  as  a 
function  of  laser  irradiance. 

II.  RESULTS  FROM  EXPERIEMTS 

In  previous  experiments,  a  Nd:YAG  laser  with  266,  532, 
and  1064  nm  wavelengths  and  3  ns  pulse  duration  is  focused 
to  —35  fjtm  diameter  spot  on  a  single-crystal  silicon  target; 
the  laser  irradiance  varies  from  3  X 109  to  1.2 
X  1012  W/cm2.  A  second  Nd-YAG  laser  of  pulse  duration  35 
ps  was  used  as  a  probe  beam  to  record  shadowgraph  images 
of  mass  ejection  from  the  laser-irradiated  silicon  surface. 
This  532  nm  wavelength  probe  beam  was  aligned  parallel  to 
the  silicon  surface  and  perpendicular  to  the  ablation  laser 
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FIG.  1.  Comparison  of  measured  ablation  depths  with  the  computational 
data  (■:  measured  ablation  depth  for  266  nm  wavelength,  0:  measured 
ablation  depth  for  532  nm  wavelength,  +:  measured  ablation  depth  for  1064 
wavelength.  The  solid  line  for  the  computational  ablation  depth  for  266  nm 
wavelength,  the  dash  line  for  the  computational  ablation  depth  for  532  nm 
wavelength,  and  the  dot  line  for  computational  ablation  depth  for  1064  nm 
wavelength.). 


beam.  The  time  delay  of  the  probe  laser  pulse  with  respect  to 
the  ablation  pulse  was  controlled  by  a  four-channel  digital 
delay/pulse  generator  and  monitored  using  a  digital  oscillo¬ 
scope.  The  depth  of  the  crater  was  measured  with  a  white 
light  interferometric  microscope  (Zygo  model  NewView 
200).  By  measuring  ablation  depth  versus  laser  irradiance  at 
different  wavelengths  (Fig.  1),  it  was  realized  that  the  abla¬ 
tion  depth  increases  dramatically  at  the  laser  irradiance 
threshold  of  about  2  X  1010  and  4  X  1011  W/cm2  for  266  and 
532  nm  wavelength  laser  pulses,  respectively.  For  the  266 
nm  wavelength,  the  ablation  depth  increases  gradually  from 
0.6  to  1.5  pm  as  the  laser  irradiance  increases  from  3.0 
X  109  to  2.0 X  1010  W/cm2,  then  it  abruptly  increases  from 
1.5  to  6.3  pm,  and  reaches  20  pm  at  1.5  X1011  W/cm2. 
For  532  nm  wavelength,  the  ablation  depth  increases  from 
0.26  to  2.5  pm  as  laser  irradiance  is  increased  from  3.4 
X  109  to  4.0  X  1011  W/cm2,  and  beyond  a  threshold,  the  ab¬ 
lation  depth  changes  sharply  from  2.5  to  7.4  /im.  For  1064 
nm  wavelength  laser  pulse,  however,  such  a  dramatic  in¬ 
crease  does  not  exist,  the  ablation  depth  increases  gradually 
from  0.24  to  3.8  /rm  when  laser  irradiance  varied  from  3.4 
X  109  to  5.6X1011  W/cm2,  then  gradually  decreases  to 
1.7  pm  at  1.2  X  1012  W/cm2.  A  shock  wave  was  observed 
to  lasts  about  several  tens  of  nanoseconds  after  the  laser 
pulse,  which  was  formed  due  to  the  pressure  difference  be- 
tween  a  dense  plasma  and  the  ambient.  For  266  and  532 
nm  wavelength  laser  pulses,  when  the  laser  irradiance  ex¬ 
ceeds  the  respective  threshold,  there  are  large  particulates 
ejected  from  the  target  surface  about  300-400  ns  after  the 
laser  irradiation,  while  for  1064  nm  wavelength,  there  were 
no  such  large  particulates. 


III.  ANALYSIS 

There  are  two  basic  types  of  thermal  processes  for  laser 
ablation/  The  first  is  normal  vaporization,  which  occurs  at 
extreme  outer  surface  for  any  laser  irradiance  and  pulse 
duration.14'15  The  vaporization  flux  is  governed  by  the  Hertz- 


Knudsen  equations,  and  the  velocity  of  surface  recession  can 
be  given  if  multiplied  by  ml p  (m  being  the  particle  mass  and 
p  the  target  mass  density): 


dx 

jf=o 


/3(Pb  ~  Famb)  ~(2mnkBT)  1/2 
P 


Xexp 


Levm 

1 

1)1 

[  kB 

\Tb 

77  J 

cm/s. 


(1) 


Here,  fi  is  the  vaporization  coefficient,16  pb  is  the  boiling 
pressure  (—0.1  MPa),  pamb  is  the  ambient  vapor  pressure,  Th 
is  the  corresponding  boiling  temperature,  kB  is  the  Boltz¬ 
mann  constant,  and  Lcy  is  the  heat  of  vaporization. 

The  second  mechanism  for  laser  ablation  is  explosive 
boiling,  much  of  the  theoretical  foundation  on  explosive 
boiling  was  established  by  Martynyuk.  ’  The  theory  of  ex¬ 
plosive  boiling  may  be  considered  from  either  a  thermody¬ 
namic  or  kinetic  viewpoint.19"0  The  former  provides  a  rigor¬ 
ous  method  to  predict  the  thermodynamic  critical 
temperature,  while  the  latter  mechanism  models  the  rate  of 
formation  of  vapor  bubble  growth  at  any  temperature.  Ac¬ 
cording  to  thermodynamic  theory  of  explosive  boiling,  the 
liquid  begins  to  be  superheated  and  becomes  metastable 
when  it  exceeds  a  temperature  of  about  ().807'tc.  When  the 
laser  irradiance  is  high  enough  (109  W/cm2  or  above),  the 
temperature  of  the  target  heated  by  the  laser  pulse  can  exceed 
0.80 Ttc,  as  a  result,  a  superheated  metastable  liquid  layer 
forms.  Homogeneous  bubble  nucleation  will  occur  in  this 
superheated  layer,  and  the  “liquid”  is  essentially  a  mixture  of 
liquid  droplets  and  vapor  that  can  facilitate  explosive  boiling. 

It  is  reasonable  to  suggest  that  mass  below  the  irradiance 
threshold  is  removed  at  the  surface  by  normal  vaporization, 
while  dramatic  increase  in  ablation  depth  is  related  to  the 
explosive  boiling.  However,  whether  explosive  boiling  oc¬ 
curs  during  the  laser  pulse  or  after  the  laser  pulse  has  not 
been  fully  established?  Miotello  and  Kelly  argued  that  ex¬ 
plosive  boiling  may  occur  during  the  laser  pulse.  However, 
in  our  case,  from  the  experiment  as  well  as  kinetic  theory  of 
explosive  boiling,  we  found  that  explosive  boiling  should 
take  place  after  the  completion  of  laser-material  interaction. 
First,  the  experimental  results  show  that  large  size  particu¬ 
lates  are  ejected  about  300-400  ns  after  laser  pulse.  Second, 
according  to  kinetic  theory  of  explosive  boiling,  it  is  not  an 
inevitable  process  when  the  liquid  is  superheated.  Homo¬ 
geneous  bubble  nucleation  occurs  when  the  liquid  is  super¬ 
heated,  but  only  if  these  bubbles  have  enough  time  to  reach  a 
critical  radius  rc,  will  they  grow  spontaneously.  Then  the 
liquid  experiences  large  density  fluctuation,  and  large  size 
particulates  are  then  ejected  from  the  target  surface.  The 
critical  radius  rc  can  be  expressed  as22 


2cr 

P^t(Ti)e\p{v,[p,  -  psat(Ti)]/RvT,}  -  p, ' 


(2) 


where  cr  is  the  surface  tension  of  the  liquid,  Rv  is  the  gas 
constant,  and  p;  is  the  densities  of  superheated  liquid/vapor 
with  Vf  =  1 1  pi-  T,  is  the  temperature  of  the  superheated  liquid, 
which  can  be  taken  as  0.8 5Ttc  when  explosive  boiling  occurs. 
Using  the  method  suggested  by  Martynyuk,  we  estimated 
that  the  thermodynamic  critical  temperature  of  silicon  is  ap- 
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proximately  5200  K.  psM  is  the  saturation  pressure  at  the 
superheated  liquid  temperature,  which  can  be  obtained  from 
the  Clausius-Clayperon  relation.  pt  is  the  pressure  of  the 
superheated  liquid,  and  it  can  be  approximated  by  the  recoil 
pressure  of  the  evaporating  vapor,  which  is  0.54/?sat(7'/).23 
According  to  the  power  law  relation  of  surface  tension  cr  for 
liquid  metal,  the  surface  tension  drops  about  80%  at  the  as- 
sumed  7j.  This  relation  is  also  satisfied  for  silicon  as  liquid 
silicon  behaves  much  like  a  liquid  metal."4  Using  these  pa¬ 
rameters,  we  estimate  rc  to  be  approximately  0.6  fim.  How¬ 
ever,  the  thermal  penetration  depth  during  a  laser  pulse  of 
duration  ris  xth=0.969[(kr)1/2],"5  ( k  is  the  thermal  diffusivity 
of  the  liquid  silicon,  which  is  about  0.75  cm2/s),  which  in 
our  case  is  about  0.47  /rm  for  3  ns  laser  pulse.  The  critical 
diameter  of  the  bubble  is  dc=2rc,  or  1.2  jam,  which  is  larger 
than  the  thermal  penetration  depth;  the  bubble  cannot  grow 
to  its  critical  radius  during  the  laser  pulse.  Experimental  evi¬ 
dence  suggests  that  explosive  boiling  occurs  only  if  the  su¬ 
perheated  layer  is  thick  enough.19  Even  if  the  explosive  boil¬ 
ing  occurs,  the  ablation  depth  is  approximately  equal  to  the 
thermal  penetration  depth,  which  is  about  0.47  jam.  How¬ 
ever,  in  our  experiment  the  maximum  ablation  depth  is  about 
22  /am.  Therefore,  explosive  boiling  should  occur  after  the 
laser  pulse  when  the  superheated  layer  is  sufficiently  thick. 

In  the  superheated  liquid,  bubble  growth  is  strongly  in¬ 
fluenced  by  the  momentum  interaction  between  the  growing 
bubble  and  the  surrounding  liquid  being  pushing  away.  The 
bubble  growth  in  this  time  regime  is  given  by“" 


R(t)  = 


Tl  ~  EsatW 

^evPu  1 

.  ^sat(/V)  . 

pi  J 

(3) 


pv  is  the  densities  of  superheated  liquid/vapor,  Tsat  is  satura¬ 
tion  temperature  of  the  superheated  liquid  pressure,  which 
can  be  obtained  from  the  Clausius-Clayperon  relation.  tc  is 
the  time  for  vapor  bubble  to  grow  from  initiation  (of  very 
small  radii,  several  mean  atomic  spacings)  to  the  critical  ra¬ 
dius  rc.  It  can  be  calculated  by  letting  left  side  of  Eq.  (3) 
equal  rc,  and  we  can  calculate  that  tc  is  about  70  ns.  It 
indicates  that  the  bubble  will  take  about  70  ns  to  grow  to  the 
critical  radius  of  0.6  /am,  then  the  superheated  liquid  trans¬ 
forms  into  a  mixture  of  vapor  and  liquid  droplets,  followed 
by  explosive  boiling.  However,  our  laser  pulse  duration  is 
only  3  ns,  from  Eq.  (3),  the  bubble  radius  can  only  grow  up 
to  about  0.025  /rm  in  3  ns.  As  a  result,  without  efficient 
energy  dissipation,  the  liquid  temperature  can  exceed  the 
critical  temperature  if  the  laser  irradiance  is  sufficiently  high. 
Therefore,  it  is  unlikely  that  explosive  boiling  will  occur 
during  laser  pulse.  A  high  temperature  region  will  be  formed 
near  the  surface,  which  will  eventually  penetrates  into  the 
target.  The  bubbles  will  reach  the  critical  radius  after  a  time 
of  tc;  in  experiments,  micron-sized  droplet  ejection  occurred 
300-400  ns  after  the  completion  of  the  laser  pulse.  Violent 
ejection  of  particulates  lasts  for  several  microseconds;  For 
this  period  of  time,  we  can  estimate  the  thermal  penetration 
depth  is  on  the  order  of  10  jam,  which  is  consistent  with  the 
experimental  ablation  depth  for  laser  irradiance  above  the 
threshold. 
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FIG.  2.  The  processes  of  laser  ablation — explosive  boiling. 


From  the  kinetic  theory,  we  can  also  calculate  the  rate  of 
homogeneous  nucleation  by-6 

/„  ~  1.5  X  1032  exp(-  AG„/kBT)exp(-  rhn/f)  nuclei/cm3  s 

(4) 

Here  AG„  is  the  free  energy  for  the  formation  of  a  stable 
homogeneous  nucleus.  Thrl  is  the  relevant  time  constant, 

27 

which  ranges  from  1  to  100  ns  according  to  Martynyuk, 
and  as  the  estimation  of  Kelly  and  Miotello  it  is  about  30 
ns.  /„  is  numerically  significant  (i.e.,  /„  >  I )  only  when  the 
temperature  is  near  to  Ttc.  As  an  example,  the  value  for  Cs  is 
/„=1  nucleus / cm3  s  at  7=0.8747'tc  and  7„ 
=  1026  nuclei/cm3  s  at  7  =0.9057’tc.  The  number  of  homoge¬ 
neous  nuclei  which  would  be  generated  during  the  laser 
pulse  is  7„V*r,  where  V*  is  the  heated  volume  during  the 
laser  pulse,  and  V*=^i rxthD2aser.  Z)laser  is  the  width  of  the 
laser  pulse  (in  our  experiments,  35  /am),  so  V*~4.5 
X  10-10  cm3.  Take  7„=1026  nuclei/cm3  and  rhll  =  50  ns,  the 
homogeneous  nuclei  generated  during  the  laser  pulse  equals 
5.  For  such  a  low  generation  ration  of  nuclei,  explosive  boil¬ 
ing  cannot  be  expected  to  occur  during  laser  pulse. 

From  the  above  analysis,  very  few  vapor  bubble  can  be 
generated  near  the  surface  of  the  target  during  the  laser 
pulse.  Even  if  the  bubbles  are  formed,  they  do  not  have 
enough  time  to  grow  up  to  the  critical  radius,  thus  explosive 
boiling  will  not  occur.  Therefore,  little  energy  provided  by 
the  laser  pulse  can  be  dissipated,  and  a  high  temperature 
layer  will  be  formed  at  and  below  the  target  surface  during 
the  laser  pulse  with  a  depth  approximately  equal  to  the  ther¬ 
mal  penetration  depth.  At  the  same  time  the  target  undergoes 
normal  vaporization  from  the  extreme  outer  surface,  and  an 
expanding  plasma  will  form  due  to  the  interaction  of  laser 
pulse  with  the  vapor.  As  a  consequence,  a  portion  of  the  laser 
energy  will  be  shielded  by  the  plasma.  Mass  ablation  below 
the  laser  irradiance  threshold  is  dominated  by  this  normal 
vaporization  mechanism.  The  vaporization  flux  is  governed 
by  Eq.  (1).  At  high  laser  irradiance,  after  the  laser  pulse  is 
completed,  the  high  temperature  liquid  layer  will  propagate 
into  the  target  with  thermal  diffusion.  New  bubbles  will 
emerge  in  the  superheated  liquid,  after  the  thickness  of  the 
superheated  layer  is  larger  than  the  critical  diameter,  the 
bubbles  will  grow  spontaneously  inside  it.  This  superheated 
layer  will  eventually  leave  the  target.  This  process  is  de¬ 
scribed  in  Fig.  2.  Consistent  with  experiments,  large  particu- 
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lates  are  ejected  from  the  target  several  hundred  nanoseconds 
after  laser  pulse. 


B=-nekBTe  +  nekB0i. 


(7) 


IV.  NUMERICAL  SIMULATIONS 


A  numerical  model  based  on  this  diffusion-phase  explo¬ 
sion  mechanism  has  been  established  to  determine  the  depth 
of  the  superheated  layer.  The  computed  depths  are  used  as  an 
indication  of  the  liquid  layer  that  is  indicative  to  explosive 
boiling.  To  solve  the  time-dependent  temperature  distribution 
in  the  sample,  the  following  conservation  of  energy  equation 
for  laser  heating  is  used: 


dT  d  (  dT\ 

pC~  =  —  k—  ;  +  a/laser  exp(-  ax), 
dt  ox\  ox  J 


(5) 


The  first  term  is  the  kinetic  energy  and  the  second  is  ioniza¬ 
tion  energy.  Te  and  ne  are  electron  temperature  and  density.  9, 
is  the  ionization  potential,  which  is  about  94  584  K.  The  first 
term  is  kinetic  energy  and  the  second  is  ionization  energy. 
The  energy  of  particles  is 


r.  ftpkBTp  5 


(8) 


where  Tp  and  np  are  particle  temperature  and  density,  respec¬ 
tively.  Here,  particles  include  ions  and  neutral  atoms,  so  the 
particle  density  can  be  expressed  as 


where  T  is  the  temperature,  C  is  the  specific  heat,  and  7laser  is 
the  laser  irradiance  that  reaches  the  surface  of  the  silicon 
target.  The  spatial  coordinate  x  is  in  the  direction  normal  to 
the  sample  surface  with  the  origin  located  at  the  surface. 

For  solving  Eq.  (5),  boundary  conditions  are  required  at 
a  =  0  and  x=L,  where  L  is  the  length  of  the  computational 
domain.  At  x=L,  the  temperature  of  the  material  is  assumed 
to  be  unaffected  by  the  laser  irradiation,  i.e.,  T( L,t>0)  =  70, 
T0  is  the  initial  temperature  of  the  solid.  For  x=0,  energy  loss 
to  the  surrounding  air  can  be  ignored  because  the  heat  flux  is 
negligible  compared  to  the  latent  heat  of  the  evaporating 
vapor.  Thus,  an  adiabatic  boundary  condition,  i.e.,  dTI  dx= 0, 
was  used  when  there  is  no  surface  evaporation.  When  vapor¬ 
ization  exists  at  the  surface,  the  energy  loss  was  calculated  as 
the  latent  heat  of  evaporating  mass.  The  energy  reaching  the 
surface  of  the  target  7laser,  is  then  equal  to  the  difference 
between  the  incident  laser  energy  and  the  latent  heat  of 
evaporation  at  the  surface, 

7iascr  =  (l  -R)I{Q,t)-LmPlv,  (6) 

where  v  is  the  receding  velocity  of  the  sample  surface  during 
evaporation,  R  is  the  reflectivity  of  laser  energy  at  the  target 
surface,  and  7(0,  t)  is  the  laser  irradiance  at  the  surface  of  the 
target. 

We  use  a  Gaussian  energy  distribution  as  the  laser  pulse 
source.  The  spatial  domain  of  60  or  120  /mm  (depending  on 
the  wavelength  of  the  laser)  was  uniformly  divided  with  grid 
size  of  15  nm.  The  time  step  is  1  ps.  The  initial  temperature 
of  the  sample  is  300  K.  The  thermal  and  optical  properties  of 
silicon  for  different  laser  wavelengths  have  been  published  in 
literature.  29-34 

We  include  in  the  model  the  absorption  of  laser¬ 
generated  plasma  from  the  target  surface.  Such  a  plasma  is 
frequently  observed  during  high  power  laser  ablation  of  sol¬ 
ids.  However,  it  has  previously  not  been  included  for  mod¬ 
eling  laser  ablation  in  the  explosive  boiling  regime.  The 
laser-induced  plasma  above  the  surface  can  shield  a  portion 
of  the  laser  beam  during  the  pulse,  here  we  only  consider  the 
loss  of  laser  energy  due  to  the  collision  by  electrons  with 
atoms  and  ions  (inverse  Bremsstrahlung  process).34’35  In  this 
model,  the  ablated  vapor  is  described  by  three  components: 
electrons,  ions,  and  neutral  atoms  with  subscripts  e ,  i,  and  a, 
respectively. 

35—37 

The  energy  of  electrons  is 


np  =  ni  +  na,  (9) 

where  n,  and  na  are  the  ion  and  atom  densities,  respectively. 
We  assume  that  the  electron  density  is  equal  to  the  ion  den¬ 
sity,  i.e.,  ne=iij. 

The  change  in  electron  energy  equals  energy  absorbed 
from  the  laser  pulse  minus  energy  transferred  to  particles  via 
collision.  Assuming  that  the  density  and  temperature  are  ho¬ 
mogeneous  in  the  entire  vapor  cloud,  the  energy  equation  of 
electrons  can  be  written  as 


dee 

dt 


[1  -exp(-  Hk[)] 


77 


(1  -  7?0)70(f)  -  ~kB(Te  -  Tp)vttne, 


(10) 


where  70(f)  is  the  laser  intensity  as  a  function  of  time,  vtr  is 
the  energy  transfer  frequency,  kl  is  the  absorption  coefficient 
of  the  plasma,  77  is  the  thickness  of  plasma,  and  R0  is  the 
reflectivity  of  the  target.  The  reflectivity  decreases  to  zero 
during  the  laser  pulse  after  surface  melting  start,  so  here  we 
choose  7?0=0.  The  plume  is  assumed  to  expand  with  sonic 
velocity, 

dH_  j  ykBT„ 

dt  \  M  ’ 


where  y  is  the  ratio  of  Cp  and  Cv,  and  M  is  particle  mass. 
The  absorption  coefficient  for  inverse  Bremsstrahlung  is  de¬ 
fined  as37 


k  i  = 


3.7  X  108 
a>\tre 


hu> 


- 1 


+ 


e~vene 

7T7  neat 


(ID 


where  a>  is  the  laser  frequency,  m  is  the  electron  mass,  e  is 
the  electron  charge,  and  c  is  the  light  speed.  vc  is  the 
electron-atom  collision  frequency, 


8kBTP 


(12) 


where  crc  is  the  electron-atom  collision  cross  section. 

The  change  in  particle  energy  is  equal  to  energy  carried 
by  the  atom  vapor  ejected  from  the  sample  surface  plus  the 
energy  transferred  from  electron  via  collision. 
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deD  3  Jv  3 

— f +  -fcB(re  -  T„)vtIne,  (13) 

dt  2  H  2  p 

where  T,  is  the  surface  temperature  of  the  sample. 

•  37 

The  ionization  rate  is 

dne  2 

—  =  a,nane- /3Rnlne,  (14) 

dt 

where  al  =  ClSkljTe  exp(- 0,/  Te)  and  /3R 

=  C,l2(2inn/ h)il2(gjl  gk)kBTe],  C,  is  an  experimental  param¬ 
eter  determined  from  electron  ionization  cross  section  data. 
The  details  of  the  calculations  of  the  parameters  utr,  ah  and 
/ 3r  can  be  found  in  Ref.  37. 

The  change  in  the  density  of  atoms  in  the  plasma  is 
equal  to  atoms  evaporated  from  the  surface 

d[H(na  +  n,)] 

dt  ~  v'  U  J 

where  Jv  is  the  evaporation  rate,  which  can  be  calculated  by 
multiplying  evaporation  flux  and  number  density  of  atom  in 
the  liquid. 

The  laser  irradiance  transmits  through  plasma  and  reach¬ 
ing  the  target  surface  7laser  can  be  written  as 

/(0,t)  =  /0(t)exp(— /ft,).  (16) 

The  ablation  depth  by  evaporation  was  calculated  by  in¬ 
tegrating  Eq.  (4).  During  the  laser  pulse,  a  high  temperature 
layer  is  formed  at  and  beneath  the  surface  of  the  target;  this 
layer  then  propagates  into  the  target  by  thermal  diffusion.  We 
consider  the  liquid  whose  temperature  is  larger  than  0.807’tc 
as  superheated  liquid,  and  in  such  a  metastable  state,  homo¬ 
geneous  bubble  nucleation  will  occur,  if  the  thickness  of  the 
superheated  layer  is  large  enough  (—4.0  /rm  according  to 
experiments).  Ablation  for  irradiance  below  the  threshold  is 
governed  by  normal  evaporation  according  to  Eq.  (4).  Abla¬ 
tion  for  the  irradiance  larger  than  the  threshold  is  governed 
by  both  normal  evaporation  and  explosive  boiling. 

The  ablation  depths  predicted  by  this  model  are  com¬ 
pared  to  experimental  data  in  Fig.  1.  The  theoretical  calcula¬ 
tions  for  all  three  wavelengths  are  consistent  with  experi¬ 
ments.  The  model  predicts  that  the  laser  irradiance  threshold 
for  explosive  boiling  is  about  3  X  101()  and  1.5 
X1011  W/cm2  for  266  and  532  nm  wavelengths,  respec¬ 
tively,  in  good  agreement  with  experimental  values  (—2 
X  1010  and  — 4X  1011  W/cm2  for  266  and  532  nm  wave¬ 
lengths,  respectively).  In  our  model,  the  superheated  liquid 
reaches  its  maximum  depth  of  100-1000  ns  after  the  laser 
pulse  is  completed,  which  also  agrees  with  experimental  re¬ 
sults.  For  1064  nm  wavelength  laser  ablation  of  silicon,  from 
theoretical  calculation,  explosive  boiling  also  does  not  occur. 

The  temperature  distribution  at  different  times  is  de¬ 
scribed  in  Fig.  3.  For  266  and  532  nm  wavelengths,  the  laser 
irradiances  are  chosen  at  their  thresholds  of  explosive  boil¬ 
ing,  which  are  3  X  10 10  and  1.5  X  1011  W/cm2,  respectively. 
For  1064  nm  wavelength,  the  irradiance  for  maximum  abla¬ 
tion  depth  is  used.  A  high  temperature  layer  is  formed  during 
the  laser  pulse.  After  the  laser  pulse,  the  high  temperature 
layer  penetrates  into  the  target,  and  the  surface  temperature  is 
not  the  highest  due  to  evaporation  at  the  surface.  For  266  nm 


0  3  6  9  12 

Distance(i^m) 


0  3  6  9  12 

Distance(pm) 

FIG.  3.  The  temperature  distribution  at  different  times,  (a)  266  nm  wave¬ 
length,  laser  irradiance  of  1=3  X  1()10  W/cm2.  (b)  532  nm  wavelength,  laser 
irradiance  of  7=1.5X10n  W/cm2.  (c)  1064  nm  wavelength,  laser  irradi¬ 
ance  of  7=4.0 X  1011  W/cm2.  (The  layer  above  the  solid  line  is  superheated 
layer). 


laser  pulse,  the  superheated  layer  reaches  its  maximum  thick¬ 
ness  about  4.0  /jl m  at  75  ns.  For  532  nm  wavelength,  the 
superheated  layer  has  a  maximum  thickness  of  4.5  /xm  at 
100  ns.  For  1064  nm,  it  reaches  a  maximum  superheated 
layer  at  30  ns,  with  a  thickness  of  3.2  /xm,  then  it  reduces 
gradually,  at  125  ns  the  superheated  layer  vanishes,  this  can 
explain  why  explosive  boiling  does  not  occur. 

Plasma  shielding  plays  an  important  role  in  determining 
the  laser  irradiance  threshold  for  explosive  boiling.  The  ef¬ 
fect  of  plasma  shielding  can  be  illustrated  by  plotting  the 
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Time  (ps) 


FIG.  4.  Temporal  profiles  of  laser  irradiance  on  the  target  surface  for  dif¬ 
ferent  initial  peak  laser  irradiances,  7peak,  before  the  interaction  with  a  mass 
plasma,  (a)  266  nm  wavelength,  the  values  of  7peak  are  a:  1()10,  b:  2  X  1()10,  c: 
3  X  1010,  and  d:  1  X  1011  W/cm2.  (b)  532  nm  wavelength,  the  values  of  7peak 
are  a:  1010,  b:  4X  1010,  c:  1.5  X  1011,  and  d:  4X  1011  W/cm2.  (c)  1064  nm 
wavelength,  the  values  of  7peak  are  a:  1010,  b:  5  X  1010,  c:  1  X  1011,  and  d: 
8  X  10u"w/cm2. 

transmitted  laser  temporal  profile  through  the  plasma  (Fig. 
4).  When  the  laser  irradiance  is  low,  the  laser  pulse  retains  its 
original  profile  with  little  attenuation  by  the  plasma.  How¬ 
ever,  when  the  laser  irradiance  is  large,  the  trailing  part  of  the 
laser  pulse  is  truncated.  The  larger  the  laser  irradiance,  the 
more  the  laser  energy  is  truncated. 

V.  CONCLUSIONS 

As  measured  from  experiments,  for  266  and  532  nm 
wavelength  laser  ablations  of  silicon,  there  exists  a  laser  ir¬ 


radiance  threshold  that  corresponds  to  a  dramatic  increase  in 
the  ablation  depth.  Below  the  threshold,  the  ablation  depth 
increases  slowly.  However,  at  the  threshold,  the  ablation 
depth  increases  dramatically,  which  accompanies  large  par¬ 
ticulates  ejected  300-400  ns  after  the  completion  of  ablation 
laser  pulse.  We  developed  a  theoretical  analysis  for  this 
threshold  phenomenon,  which  includes  the  effect  of  plasma 
shielding.  We  conclude  that  the  mass  is  removed  by  normal 
evaporation  for  laser  irradiance  below  the  threshold.  For  la¬ 
ser  irradiance  above  the  threshold,  during  the  laser  pulse  the 
mass  is  removed  by  normal  evaporation.  At  the  same  time,  a 
plasma  forms  due  to  the  interaction  of  laser  pulse  with  the 
vapor,  which  shields  part  of  the  laser  energy.  A  high  tempera¬ 
ture  layer  near  the  surface  of  the  target  is  formed  during  the 
laser  pulse,  which  penetrates  into  the  target  after  the  comple¬ 
tion  of  laser  interaction.  When  the  superheated  is  thick 
enough,  explosive  boiling  occurs.  Plasma  shielding  during 
laser  irradiation  was  found  to  have  a  significant  effect  on  the 
threshold  phenomenon,  and  our  calculations  provide  a  satis¬ 
factory  estimate  of  the  experimental  results.  Although  our 
experiment  and  our  model  are  based  on  3  ns  laser  pulse 
interaction  with  silicon,  it  should  be  applicable  for  a  broad 
range  of  pulse  durations  according  to  the  analysis  that  is  not 
material  specific. 
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